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This  paper  reviews  the  second  law  analysis  of  thermodynamics  in  enclosures  due  to  buoyancy-induced 
flow  for  energy  systems.  It  defines  entropy  generation  minimization  or  thermodynamic  optimization. 
In  addition,  the  article  summarizes  the  recent  works  on  entropy  generation  in  buoyancy-induced  flows 
in  cavity  and  channels.  Studies  on  mixed  convection  were  also  included  in  the  study.  Presentation  was 
performed  for  flow  in  porous  media  and  viscous  fluid  filled  media  at  different  shaped  enclosures  and 
duct  under  buoyancy-induced  force. 
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1.  Introduction 

Entropy  generation  minimization  or  thermodynamic  optimiza¬ 
tion  is  not  an  old  technique.  It  is  also  components  of  exergy  analysis 
[1].  The  method  of  thermodynamic  optimization  or  entropy  gen¬ 
eration  minimization  (EGM)  is  a  field  of  activity  at  the  interface 
between  heat  transfer,  engineering  thermodynamics,  and  fluid 
mechanics.  The  entropy  generation  minimization  technique  is 
based  on  the  simultaneous  application  of  the  first  law  and  the 
second  law  in  analysis  and  design.  Entropy  can  be  used  to  distin¬ 
guish  between  reversible  and  irreversible  processes.  It  analysis  can 
quantify  the  level  of  energy  quality.  In  this  paper  we  review  the 
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fundamentals  of  the  method,  its  current  status,  and  a  few  examples 
on  natural  convection  studies. 

The  destroyed  exergy  is  proportional  to  the  generated  entropy. 
Exergy  is  always  destroyed,  partially  or  totally:  this  is  the  state¬ 
ment  made  by  the  second  law  of  thermodynamics.  The  destroyed 
exergy,  or  the  generated  entropy  is  responsible  for  the  less-than- 
theoretical  thermodynamic  efficiency  of  the  system. 

First  law  efficiency  of  thermodynamics  in  a  heat  transfer  engi¬ 
neering  system  is  very  much  restricted.  Thus,  calculations  using  the 
second  law  of  thermodynamics,  which  is  related  to  entropy  gener¬ 
ation,  are  more  reliable  than  first  law  based  calculations.  In  almost 
all  thermal  systems,  second  law-based  efficiency  can  be  defined  in 
terms  of  the  ratio  of  actual  thermal  efficiency  under  the  same  condi¬ 
tions.  Therefore,  the  second  law  of  thermodynamics  can  be  applied 
to  investigate  the  irreversibility  in  terms  of  the  entropy  genera¬ 
tion  rate.  Determination  of  entropy  generation  is  also  important  to 
enhance  system  performance,  because  entropy  generation  is  the 
measure  of  destruction  of  available  work  of  the  system  [2-4]. 
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Nomenclature 

Be 

Bejan  number 

C 

concentration 

FTI 

fluid  friction  irreversibility 

HTI 

heat  transfer  irreversibility 

k 

thermal  conductivity 

Ns 

total  entropy  generation 

S 

entropy 

T 

temperature 

U,  V ,  w 

velocities 

x,y,z 

coordinates 

volumetric  rate  of  entropy  generation 

Greek  letters 

li 

dynamic  viscosity 

0 

irreversibility  ratio 

The  method  relies  on  the  simultaneous  application  of  principles 
of  heat  and  mass  transfer,  fluid  mechanics,  and  engineering  ther¬ 
modynamics,  in  the  pursuit  of  realistic  models  of  processes,  devices, 
and  installations. 

2.  A  brief  theory 

It  is  well  known  from  the  thermodynamics  that  heat  transfer 
in  the  system  results  in  thermodynamic  irreversibility  and  it  gen¬ 
erates  the  entropy.  The  minimization  of  entropy  is  a  method  of 
optimization  of  the  thermodynamic  imperfections  and  fluid  flow 
irreversibilities  [5].  Thus,  investigation  of  entropy  generation  and 
minimization  in  the  heating  systems  is  important  to  include  irre¬ 
versibilities  into  the  thermodynamic  analysis.  Entropy  is  a  measure 
of  chaos  and  if  there  is  entropy  in  a  system,  the  quality  of  energy 
decreases.  A  system  of  high  entropy  is  more  chaotic  or  disordered 
than  one  of  low  entropy.  For  example,  a  field  with  papers  scattered 
about  has  higher  entropy  than  the  field  with  the  papers  neatly  piled 
[1  ].  A  study  of  entropy  generation  in  fundamental  convective  heat 
transfer  was  obtained  by  Bejan  [6,7]. 

Entropy  analysis  involves  the  evaluation  of  the  inputs,  outputs, 
accumulations  and  creations  of  entropy  for  a  considered  system.  It 
means  that 


By  adding  mass  transfer  the  equations  [10], 


In  above  equation,  C  defines  the  concentration.  Eq.  (3)  clearly 
shows  that  entropy  generation  in  a  convection  heat  transfer 
occur  mainly  heat  transfer  and  fluid  friction.  Thus,  an  important 
parameter  can  be  defined  as  irreversibility  ratio.  The  definition  of 
irreversibility  ratio  may  enhance  the  understanding  of  the  irre¬ 
versibilities  associated  with  the  heat  transfer  and  the  fluid  friction 
[5].  It  is  given  as  [6,7], 


S'" 

0u,fluidfriction 

=  ^ - •  (5) 

S'" 

o,  heattransfer 

Carrington  and  Sun  [38]  investigated  the  second  law  analysis 
of  combined  heat  and  mass  transfer  for  both  internal  and  exter¬ 
nal  flows.  They  found  an  equation  for  the  entropy  generation.  The 
dimensionless  Bejan  number  (Be)  [  1 1  ]  is  used  in  entropy  studies.  It 
can  be  used  as  the  alternative  irreversibility  distribution  [11], 


Be  = 


(VT)2 

(VI)2  +(V<^)2 


Be  »  0.5  is  the  limit  at  which,  heat  transfer  irreversibility  dom¬ 
inates;  Bee 0.5  is  the  opposite  limit  at  which,  irreversibility  is 
dominated  by  fluid  friction  effects;  and  Be  =  0.5  is  the  case  wherein, 
heat  transfer  irreversibility  and  fluid  flow  irreversibility  are  of  equal 
importance.  Here  also,  integrating  Bejan  numberyields  global  Bejan 
number  [12], 


^global  — 


Be(X,  Y)dXdY 


v 

Entropy  input  +  Entropy  generation  -  Entropy  output  =  Entropy  accumulation  (change  in  the  total  entropy  of  the  system)(l ) 

Other  fundamental  studies  on  entropy  generation  in  convec- 


In  above  formulation,  the  entropy  creation  is  identical  to  negen- 
tropy  consumption  and  proportional  to  exergy  consumption  [1  ]. 
Entropy  balance  in  the  rate  form  is  given  by 

ASsystem  =  (Sjnlet  —  ^outlet)  +  0  (2) 

The  local  volumetric  rate  of  entropy  generation,  S'"  (W/m3I<), 
based  on  the  second  law  of  thermodynamics,  is  given  for  three 
dimensional  flow  as  follows  [8,9], 


tive  heat  transfer  and  isothermal  convective  mass  transfer  [13], 
combined  heat  and  mass  transfer  [14],  counter  flow  heat  exchang¬ 
ers  for  gas-to-gas  applications  [15],  fundamental  convective  heat 
transfer  [16],  design  optimization  of  thermal  system  [17]  and  irre¬ 
versibilities  in  various  duct  [18].  Table  1  illustrates  a  summarize  on 
fundamental  studies  of  entropy  generation. 

3.  Studies  for  natural  convection  in  porous  media  filled 
systems 


The  first  term  on  the  right-hand  side  is  due  to  heat  transfer 
and  the  second  term  is  due  to  viscous  dissipation.  It  is  well  known 
that  the  viscous  dissipation  in  the  first  law  of  thermodynamics  is 
negligible  in  many  engineering  applications. 


Natural  convection  in  porous  media  has  an  important  topic 
due  to  the  day-by-day  improvements  of  engineering  in  industry. 
There  are  many  applications  for  this  topic  such  as  building  materi¬ 
als,  insulation  technology,  bioengineering,  geothermal  application, 
chemical  technologies,  filter  design,  etc.  [  1 9,20].  Entropy  analysis  is 
an  effective  method  for  studies  on  porous  media  to  obtain  system 
efficiency.  Fundamentals  of  entropy  generation  in  natural  convec¬ 
tion  are  defined  widely  by  Baytas  and  Baytas  [21].  Mahmud  and 
Fraser  [22]  investigated  the  nature  of  heat  transfer  and  entropy 
generation  for  natural  convection  in  a  two-dimensional  circular 
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Table  1 

Summary  of  fundamental  studies  on  entropy  generation  in  convective  heat  transfer. 


913 


Author 

System 

Topic 

Rosen  [1  ] 

Different  engineering  systems 

Second-law  analysis 

Bejan  [2] 

Different  engineering  systems 

Second-law  analysis  in  heat  transfer 

Bejan  [6-9] 

Different  engineering  systems 

Fundamentals  of  exergy  analysis,  entropy  generation  minimization 

Arpaci  and  Selamet  [27] 

Energy  systems 

Entropy  efficiency 

Carrington  and  Sun  [38] 

Combined  heat  and  mass  transfer  in  both  internal  and 
external  flows 

Second-law  analysis 

Bejan  [15] 

Heat  exchanger  design  counter  flow  heat  exchangers  for 
gas-to-gas  applications 

Irreversibility 

Tsatsaronis  et  al.  [17] 

Design  optimization 

Exergy 

Naterer  and  Camberos 

Simulation 

Entropy  and  the  second  law 

Prigogine  [63] 

Thermodynamics  processes 

Irreversibilities 

Oztop  [64] 

Semicircular  duct 

Second  law  analysis 

Amani  and  Nobari  [73] 

Curved  pipes 

Entropy  generation 

Chen  et  al.  [67] 

Vertical  channel 

Entropy  generation 

Hernandez  et  al.  [74] 

Fuel  cell 

Entropy  generation  analysis 

section  enclosure  filled  with  porous  media.  In  their  work,  the  Darcy 
momentum  equation  is  used  to  model  the  porous  media.  Results 
are  presented  in  terms  of  Nusselt  number,  entropy  generation  num¬ 
ber,  and  Bejan  number.  They  found  that  for  conduction  regime, 
both  average  Nusselt  number  and  the  entropy  generation  number 
are  independent  of  Rayleigh  number  variation.  In  convection- 
dominated  regime,  these  parameters  show  an  increasing  tendency 
with  increasing  Rayleigh  number.  At  high  Rayleigh  number,  the 
near-wall  magnitude  of  overall  entropy  generation  rate  is  higher, 
but  heat  transfer  irreversibility  is  higher  at  the  center  portion  of  the 
cavity.  These  are  summarized  in  Fig.  1.  By  plotting  Nusselt  number 
and  entropy  generation  number  versus  Rayleigh  number.  In  their 
another  study,  the  problem  of  entropy  generation  in  a  fluid  sat¬ 
urated  porous  cavity  for  laminar  magnetohydrodynamic  natural 
convection  heat  transfer  is  analyzed  by  considering  Darcy’s  law  for 
porous  media.  Magnetic  force  is  assumed  acting  along  the  direc¬ 
tion  of  the  gravity  force.  As  boundary  conditions  of  the  cavity,  two 
vertical  opposite  walls  are  kept  at  constant  but  different  temper¬ 
atures  and  the  remaining  two  walls  are  kept  thermally  insulated. 
For  a  range  of  Rayleigh  number  (Ra  =  1-104)  and  Hartmann  num¬ 
ber  (Ha  =  0-10),  heat  transfer,  overall  entropy  generation  rate,  and 
heat  transfer  irreversibility  are  presented  in  terms  of  dimension¬ 
less  Nusselt  number  (Nu),  entropy  generation  number  (Ns),  and 
Bejan  number  (Be),  respectively.  Fig.  2,  in  their  paper,  presents  Nsav 
and  Beav  for  Fig.  2(a)  and  (b),  respectively.  These  are  function  of 


Fig.  1.  Average  Nusselt  and  entropy  generation  numbers  as  a  function  of  Rayleigh 
number  [5]. 


Hartmann  number  at  different  Ra.  Similar  to  the  Nuav-Ha  profile 
characteristics,  Nsav  is  maximum  at  Ha  =  0,  decreases  with  increas¬ 
ing  Ha  and  approaches  an  asymptote  at  higher  values  of  Hartmann 
number.  Magnitudes  of  Beav  are  almost  same  for  all  Rayleigh  num¬ 
bers  considered  here  at  Ha  =  0.  Increases  in  the  value  of  Ha  have  a 
tendency  to  slowdown  the  fluid  movement  inside  the  cavity,  thus 
causing  relative  increases  of  heat  transfer  irreversibility  (i.e.,  Beav). 

Baytas  [11]  performed  a  study  to  investigate  the  entropy  gen¬ 
eration  due  to  buoyancy  forces  in  an  inclined  porous  cavity.  He 
observed  that  when  Rayleigh  number  decreases,  heat  transfer  irre¬ 
versibility  begins  to  dominate  the  fluid  friction  irreversibility.  The 
Bejan  number  is  rapidly  changed  between  1 50°  and  270°.  The  prob¬ 
lem  of  entropy  generation  for  unsteady  natural  convection  and 
radiation  in  a  tilted  open-ended  channel,  submitted  to  uniform 
hot  walls  temperature  and  filled  with  an  isotropic  fluid  saturated 
porous  medium  is  outlined  by  Slimi  [23].  In  his  study,  Darcy  flow 
model  was  used  and  viscous  dissipation  effects  are  taken  into 
account.  His  results  showed  that  the  global  heat  transfer  rate 
depends  slightly  on  the  viscous  dissipation  term  but  is  strongly  sub¬ 
ordinate  to  the  tilted  angle.  Moreover,  the  calculation  of  entropy 
generation  maps  is  feasible  and  gives  a  measure  of  the  degree 
of  irreversibility  of  the  convective  flow.  Recently,  entropy  gener¬ 
ation  due  to  buoyancy-induced  convection  and  conduction  in  a 
right  angle  trapezoidal  enclosure  filled  with  fluid  saturated  porous 
medium  has  been  performed  numerically  by  Varol  et  al.  [19].  In 
their  case,  the  left  vertical  solid  wall  of  the  trapezoidal  enclosure 
has  a  finite  thickness  and  conductivity.  The  governing  Darcy  and 
energy  equations  are  solved  numerically  using  a  finite  difference 
method.  Entropy  generation  is  calculated  by  using  the  obtained 
velocities  and  temperature  distributions  from  the  computer  code. 


Fig.  2.  (a)  Average  entropy  generation  number  as  a  function  of  Rayleigh  number, 
(b)  Average  Bejan  number  as  a  function  of  Rayleigh  number  [16]. 
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Results  are  presented  for  the  Bejan  number,  local  and  mean  Nusselt 
numbers,  streamlines,  isotherms,  iso-Bejan  lines  and  entropy  gen¬ 
eration  contours.  It  is  found  that  the  most  important  parameters 
on  heat  transfer  and  fluid  flow  are  thermal  conductivity  ratio  and 
dimensionless  thickness  of  the  solid  wall  of  the  enclosure.  Thus, 
these  parameters  also  generate  entropy  for  the  whole  system.  It 
is  also  observed  that  increasing  the  Rayleigh  number  decreases 
the  Bejan  number;  however,  heat  transfer  is  an  increasing  func¬ 
tion  of  Rayleigh  number.  Non-uniform  heating  effects  are  tested 
on  entropy  generation  due  to  natural  convection  in  a  isosceles  tri¬ 
angular  enclosures  at  different  positions  by  Varol  et  al.  [24].  They 
found  that  both  inclination  angles  and  Rayleigh  numbers  make 
important  effect  on  natural  convection  heat  transfer,  fluid  flow  and 
entropy  generation.  The  highest  entropy  generation  due  to  HTI  and 
FFI  and  stream  function  are  observed  at  (p  =  90°.  Multiple  cells  were 
formed  at  this  angle.  Streamlines,  isotherms  and  entropy  contours 
are  symmetric  inside  the  enclosure  for  both  cp  =  0°  and  cp  =  1 80°. 

In  Zahmatkesh’s  [12]  study,  a  numerical  analyze  was  per¬ 
formed  on  the  importance  of  thermal  boundary  conditions  of  the 
heated/cooled  walls  in  heat  transfer  and  entropy  generation  char¬ 
acteristics  inside  a  porous  enclosure,  heated  from  below.  Both 
the  heating  and  the  cooling  are  carried  out  uniformly  and  non- 
uniformly  and  the  results  are  compared.  He  presents  results  in 
terms  of  streamlines,  isothermal  lines,  iso-entropy  generation  lines, 
and  iso-Bejan  lines.  Additionally,  variations  of  average  Nusselt 
number,  global  entropy  generation  rate,  and  global  Bejan  number 
are  analyzed  over  a  wide  range  of  Darcy-modified  Rayleigh  num¬ 
ber  ( 1 0  <  Ra  <  1 000).  Inspection  of  the  results  indicates  that  thermal 
boundary  conditions  are  of  profound  influences  on  the  induced 
flow  as  well  as  heat  transfer  characteristics  and  possess  prominent 
consequences  on  entropy  generation  rates.  It  is  demonstrated  that, 
the  optimum  case  with  respect  to  heat  transfer  as  well  as  entropy 
generation  could  be  achieved  by  non-uniform  heating. 

In  another  study  of  Baytas  [25],  the  entropy  generation  in  a  sat¬ 
urated  porous  cavity  for  laminar  natural  convection  is  investigated 
by  using  the  non-Darcy  flow  model  and  a  thermal  nonequilibrium 
model  of  the  heat  transfer  between  the  fluid  and  the  solid  phases. 
His  results  demonstrated  that  the  entropy  generation  due  to  vis¬ 
cous  drag  is  very  important  for  the  non-Darcy  flow  model,  and  the 
high  temperature  difference  between  solid  and  fluid  phases  shows 
the  importance  of  the  thermal  nonequilibrium  model  for  porous 
media.  He  suggested  that  entropy  generation  for  the  thermal 
nonequilibrium  and  non-Darcy  flow  models  should  be  considered 
for  engineering  design  and  thermodynamic  optimization.  Mahmud 
and  Fraser  [26]  investigated  the  nature  of  entropy  generation  for 
natural  convection  in  a  square  enclosure  filled  with  porous  media 
vibrating  sinusoidally  perpendicular  to  the  applied  temperature 
gradient  in  a  zero-gravity  field.  They  indicated  that  Gravity  oscil¬ 
lation  introduces  a  true  periodic  behavior  to  the  Nusselt  number, 
Bejan  number,  and  entropy  generation  rate.  The  periodic  response 
of  these  three  parameters  is  synchronized  with  the  forced  acceler¬ 
ation,  namely,  having  the  same  period  as  the  forced  acceleration. 
Table  2  gives  a  summary  of  the  various  theoretical  investigations 
for  entropy  generation  in  porous  media  filled  systems. 

4.  Studies  for  natural  convection  viscous  fluid  filled 
medium 

In  industrial  application,  natural  fluids  as  air  or  water  are  used 
mostly.  Thus,  in  the  analyses,  thermal  properties  of  these  fluids 
are  taken  into  account  to  the  analysis.  Arpaci  and  Selamet  [27] 
investigated  the  internal  energy  generation  in  case  of  buoyant  tur¬ 
bulent  flows.  They  developed  a  thermal  model  based  on  microscales 
flows  for  different  parameters.  Dagtekin  et  al.  [28]  performed  a 
numerical  study  to  investigate  the  second  law  analysis  due  to 


Fig.  3.  Evolution  of  the  maximum  of  the  dimensionless  total  entropy  genera¬ 
tion  versus  Rayleigh  number  and  distribution  irreversibility  ratio  (p  =  10-4,  10-3, 
3.3  x  10-3  and  10“2  [30]. 

buoyancy-induced  flow  in  a  T-shaped  enclosure.  They  investigated 
effects  of  geometrical  ratio  and  Rayleigh  number  on  entropy  gen¬ 
eration. 

Influence  of  thermodiffusion  on  the  birth  of  free  convection 
in  Rayleigh-Benard  configuration  was  investigated  by  Platten  and 
Chavepeyer  [29].  They  summarize  the  influence  of  thermodiffusion 
on  free  convection  together  with  the  way  to  use  this  mutual  influ¬ 
ence  in  order  to  experimentally  deduce  the  value  of  the  Soret  effect. 
Magherbi  et  al.  [30]  studied  the  time-dependent  natural  convec¬ 
tion  and  second  law  analysis  of  thermodynamics  in  a  differentially 
heated  enclosure  using  finite  volume  technique.  They  indicated 
that  the  total  entropy  generation  has  a  maximum  value  at  the 
onset  of  the  transient  state,  which  increases  with  the  Rayleigh  num¬ 
ber  and  the  irreversibility  distribution  ratio.  They  also  found  that 
entropy  generation  asymptotically  tends  towards  a  constant  value 
at  low  Rayleigh  numbers,  whereas  an  oscillation  of  the  entropy  gen¬ 
eration  was  observed  for  higher  Rayleigh  numbers,  before  reaching 
the  steady  state.  Graphs  of  the  maximum  of  the  total  entropy 
generation  Max  (Sr?a)  as  function  of  Rayleigh  number  at  differ¬ 
ent  irreversibility  distribution  ratio  ip  are  shown  in  Fig.  3.  At  low 
irreversibility  distribution  ratios  (as  low  as  ip<  10-4)  Max  (ST?a) 
takes  on  small  values  even  at  high  Rayleigh  numbers.  However, 
for  tp  <  10-3,  Max  (5T?a)  increases  rapidly.  Similar  observations  can 
be  made  about  the  evolution  of  the  total  entropy  generation  (ST?a) 
as  a  function  of  the  Rayleigh  number  during  the  steady  state  as  seen 
in  Fig.  4. 

Entropy  generation  in  rectangular  cavities  with  the  same  area 
but  different  aspect  ratios  is  numerically  investigated  by  Ilis  et  al. 
[31  ].  They  used  finite  difference  method  and  the  energy  and  vortic- 
ity  equations  are  solved  line  by  line  by  employing  the  ADI  method, 
whereas  the  stream  function  equation  is  solved  point  by  point.  It 
is  found  that  for  a  cavity  with  high  value  of  Rayleigh  number  (i.e., 
Ra  =  1 05 ),  the  total  entropy  generation  due  to  fluid  friction  and  total 
entropy  generation  number  increase  with  increasing  aspect  ratio, 
attain  a  maximum  and  then  decrease.  For  all  cavities  with  different 
aspect  ratio,  the  total  entropy  generation  increases  with  increasing 
Rayleigh  number,  though  the  rate  of  increase  is  different  as  seen 
in  Fig.  5.  For  Ra  =  100,  the  total  entropy  generation  increases  from 
the  minimum  aspect  ratio  (AR=1)  to  the  maximum  aspect  ratio 
(AR  =16)  with  the  same  order.  For  Ra  =  1 05  and  0  =  1 0-4,  the  total 
entropy  generation  for  the  AR  =  1 6  cavity  is  less  than  that  for  AR  =  9 
due  to  the  peak  point  of  the  fluid  friction  total  entropy  generation. 

Mourad  et  al.  [10]  report  a  numerical  determination  of  the 
entropy  generation  in  doubly  diffusive  convection  on  2D  approx¬ 
imation  in  a  square  inclined  cavity,  filled  with  a  fluid  (assumed 
to  be  a  perfect  gas  mixture).  Their  results  showed  that  the  total 
entropy  generation  increases  with  the  thermal  Grashof  number 
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Table  2 

Summary  of  investigations  in  porous  media  filled  enclosures. 
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Author 

Study 

Method 

Results  and  remarks 

Mahmud  and  Fraser  [16] 

Magnetohydrodynamic  free 
convection  in  a  porous 
enclosure 

Control  volume  based 
finite-volume  method 

In  the  absence  of  magnetic  force,  entropy  generation  rate 
is  relatively  higher  in  magnitude  near  two  vertical  walls 

Mahmud  and  Fraser  [22] 

Free  convection  inside  a 
two-dimensional  circular 
porous  enclosure 

Control  volume  based 
finite-volume  method  was 
used.  A  non-staggered  and 
non-uniform  grid  system  is 
used  with  a  higher  mesh 
density  near  the  walls 

At  high  Rayleigh  number,  the  near-wall  magnitude  of 
overall  entropy  generation  rate  is  higher,  but  heat  transfer 
irreversibility  is  higher  at  the  center  portion  of  the  cavity 

Zahmatkesh  [12] 

Free  convection  inside  a 
circular  porous  enclosure  with 
different  boundary  conditions 

Control  volume  method  of 
Patankar  [62] 

Entropy  generation  rate  is  likely  to  be  the  highest  for 
uniform  heating/cooling  and  the  lowest  for  non-uniform 
heating 

Baytas  [25] 

Porous  square  enclosure  filled 
with  a  heat-generating  solid 
phase 

Non-equilibrium  porous  model 
Finite  volume  finite  difference 
method 

The  entropy  generation  due  to  viscous  drag  is  very 
important  for  non-Darcy  flow  model  as  well  as  high 
temperature  difference 

Baytas  [11] 

Tilted  saturated  porous  cavity 

Finite  difference  control 
method  of  Patankar  [62] 

Inclination  angle  of  cavity  makes  important  effect  on 
degree  of  irreversibility 

Varol  et  al.  [19] 

Slimi  [23] 

Right-angle  trapezoidal 
enclosure 

Tilted  saturated  porous 
channel 

Finite  difference  method 

Inclination  angle  of  the  inclined  wall  of  trapezoidal 
geometry  affects  the  entropy  generation 

Entropy  generation  maps  is  feasible  and  gives  a  measure  of 
the  degree  of  irreversibility  of  the  convective  flow 

Mahmud  and  Fraser  [49] 

Square  porous  enclosure  with 
vibration 

Control  volume  based  finite 
volume  method 

At  the  lower  extreme  of  the  gravity  oscillation  NuAV  and 
iVsav  is  minimum  and  Beav  is  maximum 

and  the  buoyancy  ratio  for  moderate  Lewis  numbers.  Locally,  the 
irreversibility  due  to  heat  and  mass  transfer  are  nearly  identical  and 
are  localized  in  the  bottom  and  the  top  of  the  heated  and  the  cooled 
walls  respectively.  Entropy  generation  due  to  natural  convection 
has  been  calculated  for  three  radii  and  a  wide  range  of  Rayleigh 
numbers  for  an  isothermal  cylinder  by  Abu-Hijleh  et  al.  [32].  They 
indicated  that,  for  air,  entropy  generation  was  predominantly  due 
to  conduction  while  the  contribution  of  viscous  dissipation  was 
negligible.  As  the  Rayleigh  number  increased,  total  entropy  genera¬ 
tion  increased.  As  the  radius  of  the  cylinder  increased,  total  entropy 
generation  decreased. 


Ra 


*-AR=1 
*-AR=4 
AR=9 
-35—  AR=16 


AR=1 

AR=4 

AR=9 

AR=16 


Ra 


Fig.  4.  Variation  of  the  dimensionless  total  entropy  generation  in  steady  state  versus 
Rayleigh  number  and  distribution  irreversibility  ratio  (p  =  10-4, 10-3, 3.3  x  10-3,  and 
10“2  [30]. 


The  optimization  in  an  inclined  square  enclosure  for  minimum 
entropy  generation  was  analyzed  by  Baytas  [33].  Based  on  that 
study,  the  local  heat  transfer  irreversibility  and  the  local  fluid 
friction  irreversibility  change  by  the  inclination  angle  and  the  min¬ 
imum  entropy  generation  depends  considerably  on  the  inclination. 
Numerical  prediction  of  local  and  total  entropy  generation  rates  for 
natural  convection  of  air  in  a  vertical  channel  symmetrically  heated 
with  a  uniform  heat  flux  was  studied  by  Andreozzi  et  al.  [34].  They 
found  that  global  entropy  generation  increases  with  both  aspect 
ratio  and  Rayleigh  number  increase.  These  values  have  been  corre¬ 
lated  as  function  of  Ra  and  L/b;  the  correlation  obtained  by  means 
of  the  least  square  method  is 

log  S*  y  =  -13.2  +  1.91  logRa  +  0.045  (f)  (8) 

In  above  equations, 

r2  =  0.999  and  5  <  =-  <  20  x  103  <  Ra  <  106.  (9) 

b 

Erbay  et  al.  [35]  analyzed  the  thermodynamic  irreversibilities 
due  to  natural  convection  in  a  square  enclosure.  Erbay  et  al.  [36,37] 
investigated  entropy  generation  with  a  single  discrete  heat  source 
on  the  vertical  wall  of  a  cavity.  Their  work  focused  on  the  transient 
development  of  the  flow  and  temperature  fields  and  instantaneous 
entropy  generation  distributions  during  the  transient  stages. 

Shuja  et  al.  [5]  made  a  numerical  work  on  natural  convection  in 
a  square  cavity  with  a  heat  generating  body  by  including  entropy 
consideration.  They  solved  both  convection  and  conduction  (conju¬ 
gate)  mode  of  heat  transfer  to  simulate  cooling  of  micro-electronic 
equipments.  They  compared  the  effects  of  working  fluid  as  air  and 
water  and  found  that  the  flow  developed  in  the  cavity  due  to  nat¬ 
ural  convection  has  almost  the  same  mode  of  velocity  fields  for 
air  and  water.  Cooling  performance  of  heated  body  is  related  with 
location  of  the  body.  The  entropy  generation  due  to  heat  transfer  is 
relatively  lower  for  the  solid  body  located  in  left  bottom  and  right 
top  corners  of  the  cavity.  It  attains  high  values  in  the  case  of  fluid 
friction.  It  is  an  interesting  result  that  the  entropy  attains  higher 
values  in  the  case  of  air  as  compared  to  the  case  for  water.  It  means 
that  the  heated  solid  body  losses  more  heat  in  air  than  in  water 
for  the  indicated  simulation  conditions.  Table  3  illustrates  the  total 
entropy  generation  in  the  cavity  due  to  HTI  and  FFI.  In  The  case  of 
air  total  entropy  generation  due  to  heat  transfer  and  fluid  friction 
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Table  3 

Total  entropy  generated  in  the  cavity  for  various  cases  [45]. 


Working  Fluid 

Solid  position 

Entropy  generated  (J/K)  due  to 

Heat  transfer 

Fluid  friction 

Total 

Irreversibility  ratio  (0) 

Air 

Bottom 

2.379  x  10“9 

6.377  x  10~l0 

3.017  x  10“9 

0.2680 

Center 

3.228  x  10-9 

1.286  x  10-10 

3.357  x  10-9 

0.0399 

Top 

2.419  x  10“9 

2.797  x  10-10 

2.698  x  10-9 

0.1160 

Water 

Bottom 

3.238  x  10-11 

1.170  x  10~10 

1.494  x  10-10 

3.6100 

Center 

5.884  x  10-11 

2.020  x  10-11 

7.905  x  10-11 

0.3430 

Top 

2.560  x  10-11 

5.024  x  10-11 

7.584  x  10-11 

1.9600 

attains  higher  values  as  compared  to  that  for  water.  In  addition,  the 
total  entropy  generation  due  to  heat  transfer  has  maximum  value 
when  the  solid  body  is  at  the  center  of  the  cavity  while  the  total 
entropy  generation  due  to  fluid  friction  has  a  minimum  value  at 
this  location.  This  may  indicate  that  fluid  motion  generated  close 
to  the  solid  body  is  equally  distributed  in  the  vicinity  of  all  sur¬ 
faces  of  the  solid  body,  which  results  in  a  relatively  low  velocity 
flow  field.  On  the  other  hand,  the  increase  in  total  entropy  due  to 
heat  transfer  is  not  only  because  of  the  convection  effect,  but  the 


Fig.  5.  The  changes  of  total  entropy  generation  with  Rayleigh  number  for  cavities 
with  different  aspect  ratios  (a)  0  =  1O~4  and  (b)  <£  =  10~2  [31]. 


conduction  effect  as  well.  From  the  table,  the  irreversibility  ratio  is 
minimum  when  the  solid  body  is  located  at  the  center  of  the  cavity. 

Varol  et  al.  [39]  studied  the  forecasting  of  entropy  generation 
of  laminar  natural  convection  in  a  partially  cooled  square  cross- 
sectional  room  has  been  performed  using  support  vector  machines 
(SVM).  The  two-dimensional  room  was  modeled  as  floor  heating 
story  with  a  window.  Forecasting  of  entropy  generation  due  to  fluid 
friction  irreversibility  (FFI)  and  heat  transfer  irreversibility  (HTI) 
were  made  with  known  values  for  unknown  parameters  using  SVM. 
Thus,  calculation  time  was  extremely  reduced  and  values  were 
obtained  even  for  non-convergence  cases. 

Recently,  Mukhopadhyay  [40]  made  a  numerical  work  to  ana¬ 
lyze  entropy  generation  due  to  natural  convection  in  square 
enclosures  with  multiple  discrete  heat  sources.  He  examixed  the 
heaters  of  equal  length  and  strength,  the  effects  of  Rayleigh  num¬ 
ber  and  heater  position  on  flow  and  temperature  fields  and  local 
entropy  generation.  He  indicated  that  for  heaters  of  equal  length 
and  strength  ratios,  placement  of  the  heaters  close  to  the  mid¬ 
plane  of  the  cavity  is  most  undesirable  as  both  peak  temperature 
and  entropy  generation  are  maximized.  On  the  other  hand,  placing 
the  heaters  very  close  to  the  walls  or  at  a  dimensionless  spacing  of 
about  0.5  is  most  desirable  as  both  these  quantities  are  minimized. 

Yang  et  al.  [41  ]  performed  the  second  law  analysis  of  thermody¬ 
namics  on  the  laminar  film  condensation  of  pure  saturated  vapor 
flowing  in  the  direction  of  gravity  on  an  ellipsoid  with  variable  wall 
temperature.  They  found  that  entropy  generation  increases  with 
ellipticity.  Furthermore,  the  irreversibility  due  to  finite  tempera¬ 
ture  difference  heat  transfer  dominates  over  that  due  to  condensate 
film  flow  friction  and  the  local  entropy  generation  rate  decreases 
with  increasing  ellipticity  in  the  upper  half  of  ellipsoid.  Abu-Hijleh 
et  al.  [42]  made  a  numerical  prediction  of  entropy  generation  due 
to  natural  convection  from  a  horizontal  cylinder.  Fig.  5  shows  that 
an  increase  in  RaD,  i.e.,  a  higher  temperature  difference  between 
the  cylinder  and  surroundings,  results  in  an  increase  in  Sgen.  The 
main  conclusion  from  Fig.  5  is  that  the  use  of  larger  cylinders 
results  in  significantly  lower  entropy  generation  for  the  same  value 
of  RaD  due  to  a  lower  temperature  difference.  A  larger  cylinder 
has  a  lower  tangential  temperature  gradient,  which  reduces  the 
contribution  from  the  second  component  in  the  conduction  part 
of  Eq.  (8).  Increased  cylinder  size  results  in  an  increase  in  vis¬ 
cous  dissipation  in  Eq.  (8).  It  should,  however,  be  noted  that  the 
second  term  of  Eq.  (8)  was  negligible  even  at  the  largest  r0  for 
the  greatest  RaD.  Thus,  there  was  no  appreciable  increase  in  Sgen 
at  large  r0  due  to  viscous  dissipation.  Analysis  of  local  entropy- 
generation  profiles  showed  that  the  location  of  maximum  entropy 
generation  depends  on  the  size  of  the  cylinder.  For  ro  =  10~3m, 
maximum  entropy  generation  occurred  in  the  plume  region,  at 
the  top  of  the  cylinder.  For  r0  =  10-1  m,  maximum  entropy  genera¬ 
tion  occurred  in  the  inflow  region,  at  the  bottom  of  the  cylinder. 
For  10-1  >r0>  10~3,  the  contributions  from  both  the  inflow  and 
plum  regions  were  of  the  same  order.  We  may  manipulate  the  flow 
in  regions  of  low  local  entropy  generation  with  minimal  adverse 
effects  on  the  total  irreversibility  of  the  component.  Jery  et  al.  [43] 
performed  a  study  on  influence  of  an  external  oriented  magnetic 
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Table  4 

Summary  of  the  calculation  of  entropy  due  to  natural  convection  in  different  enclosures. 


917 


Authors 

Geometry 

Method 

Results  and  remarks 

Abu-Hijleh  and  Heilen  [46] 

Heated  rotating  cylinder 

Elliptic  system  of  PDEs  are 
solved  by  explicit  hybrid 
scheme 

The  total  entropy  generation  increases  with  the  increase  of 
both  Reynolds  number  and  buoyancy  parameter  but 
decreases  as  the  cylinder  radius  increases 

Abu-Hijleh  et  al.  [32] 

Heated  isothermal  cylinder 
in  oil 

Elliptic  system  of  PDEs  are 
solved  by  explicit  hybrid 
scheme 

Including  the  viscous  dissipation  term  had  minimal  effect 
on  Nusselt  number  and  total  entropy  generation 

Abu-Hijleh  et  al.  [42] 

Natural  convection  from  a 
horizontal  cylinder 

Equations  are  solved  in 

streamfunction-vorticity 

formulation 

As  the  radius  of  the  cylinder  increased,  total  entropy 
generation  decreased. 

Andreozzi  et  al.  [34] 

Vertical  channel 
symmetrically  heated  at 
uniform  heat  flux 

Numerical  solutions  of  the 
elliptic  momentum  and 
energy  equations  are 
carried  out  with  the 
stream-function-vorticity 
method 

Global  entropy  generation  increases  with  both  aspect  ratio 
and  Rayleigh  number  increase 

Baytas  [33] 

Natural  convection  in 
inclined  cavity 

Finite  volume  technique 

An  optimization  has  been  made  to  obtain  minimum 
entropy  generation 

Mukhopadhyay  [40] 

Natural  convection  in 
square  enclosures  with  2 
multiple  discrete  heat 
sources 

The  governing  equations 
were  discretized  using 
SIMPLER  algorithm  of 

Patankar  [62] 

The  dominant  contribution  to  entropy  generation  comes 
from  heat  transfer  irreversibilities,  with  fluid  friction 
accounting  for  only  a  small  fraction  of  the  total  entropy 
generation 

Famouri  and  Hooman  [52] 

Natural  convection  by 
heated  partitions  in  a 
cavity 

FORTRAN  code  and  the 
commercially  available 
CFD-ACE  software. 

While  fluid  friction  term  has  nearly  no  contribution  to 
entropy  production,  the  heat  transfer  irreversibility 
increases  monotonically  with  the  Nusselt  number  and  the 
dimensionless  temperature  difference 

ilis  et  al.  [31  ] 

Rectangular  cavity  with 
differentially  heated 
vertical  walls 

Finite  difference  technique 

For  a  cavity  with  high  value  of  Rayleigh  number  (i.e., 

Ra  =  105),  the  total  entropy  generation  due  to  fluid  friction 
and  total  entropy  generation  number  increase  with 
increasing  aspect  ratio,  attain  a  maximum  and  then 
decrease. 

Magherbi  et  al.  [30] 

Entropy  generation  at  the 
onset  of  natural  convection 
in  a  differentially  heated 
enclosure 

Of  the  control  volume 
finite-element  method 
(CVFEM) 

The  total  entropy  generation  has  a  maximum  value  at  the 
onset  of  the  transient  state,  which  increases  with  the 
Rayleigh  number  and  the  irreversibility  distribution  ratio 

Mourad  et  al.  [10] 

Numerical  determination 
of  the  entropy  generation 
in  doubly  diffusive 
convection  on  2D 
approximation  in  a  square 
inclined  cavity 

Control  volume 
finite-element  method 
(CVFEM) 

Total  entropy  generation  increases  with  the  thermal 

Grashof  number  and  the  buoyancy  ratio  for  moderate 

Lewis  numbers 

Erbay  et  al.  [37] 

Shuja  et  al.  [5] 

An  analysis  of  the  entropy 
generation  in  a  square 
enclosure  with  partial 
heating 

Square  cavity  with  heat 
generating  body 

Finite  volume  method 
(FVM)  coupled  with  power 
law  scheme  for  the 
convective  terms 

The  irreversibilities  are  dominant  due  to  heat  transfer 
whereas  fluid  friction  irreversibilities  have  been  found 
negligible  as  it  is  expected  for  the  natural  convection 

The  entropy  generation  due  to  heat  transfer  is  relatively 
lower 

Mahmud  and  Sadrul  Islam  [53] 

Wavy  enclosure 

Finite  volume  method 

Wave  ratio  is  an  important  parameter  on  entropy 
production 

Varol  et  al.  [54] 

Partially  heated  isosceles 
triangular  enclosures 

Finite  difference  method 

Length  of  partial  heater  the  most  effective  parameter  on 
entropy  generation 

Varol  et  al.  [55] 

Conjugate  natural 
convection  in  enclosures 

Finite  difference  method 

Thermal  conductivity  ratios  and  thickness  of  solid  wall 
become  effective  on  entropy  production 

Mahmud  and  Fraser  [56] 

Mixed 

convection-radiation 
interaction  in  a  vertical 
channel 

Analytical  study  with 
reasonable  simplifications 

Irreversibility  due  to  the  heat  transfer  and/or  fluid  friction 
become  zero  at  idle  points 

Bassam  [58] 

Cylinder 

Streamfunction-vorticity 

Entropy  increases  with  Ra  number 

Yilbas  [57] 

Cavity  heated  from  bottom 

Finite  volume  method 

Temperature  difference  play  important  role  on  entropy 
generation 

Mahmud  et  al.  [59] 

Vertical  channel 

Analytical  solution 

For  positive  value  of  heat  generation/absorption 
parameter,  entropy  generation  rate  is  higher  than  the 
negative  value  of  same  magnitude 

Oliveski  et  al.  [69] 

Enclosures 

General  methods 

Entropy  generation  is  applied  for  natural  convection 

Kaluri  and  Basak  [68] 

Differentially  heated 
cavities 

Finite  element  method 

Prandtl  number  plays  important  role  on  entropy 
generation 

Saleem  et  al.  [70] 

Marangoni  convection 

Finite  difference  method 

Entropy  generation  rate  increases  with  the  increase  in  the 
Marangoni  number 

Kaluri  and  Basak  [71  ] 

Discretely  heated  square 
cavity 

Finite  element  method 

High  thermal  mixing  may  not  be  the  optimal  situation  for 
achieving  uniform  temperature  distribution  based  on 
entropy  production 

Basak  et  al.  [72] 

Wide  variety  of  thermal 
boundary  conditions 

Galerkin  finite  element 
method 

Thermal  boundary  conditions  are  affected  on  thermal 
mixing,  heat  transfer  and  entropy  generation 
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field  on  entropy  generation  in  natural  convection  for  air  and  liq¬ 
uid.  They  studied  the  steady-unsteady  states  mass,  momentum  and 
entropy  equation  for  different  inclination  angle  of  the  magnetic 
field  and  Pradntl  number.  At  local  level  and  for  relatively  higher 
thermal  Grashof  number  (GrT  =  105),  entropy  generation  distribu¬ 
tion  is  strongly  dependent  on  magnetic  field  direction,  magnitude 
of  irreversibility  lines  increases  up  to  30°,  then  gradually  decreases. 
No  entropy  is  generated  in  the  cavity  center.  Table  4  summarizes 
the  various  investigations  for  entropy  generation  in  viscous  fluid 
filled  systems. 

5.  Entropy  generation  due  to  mixed  convection 

Mixed  convection  heat  transfer  occurs  for  both  natural  and 
forced  convection  present  in  a  system.  In  this  case,  Richardson 
number  (Gr/Re2)  becomes  effective.  Entropy  places  important  role 
in  this  type  of  heat  transfer  regime  due  to  heat  transfer  and  fluid 
friction. 

Narusawa  [44]  studied  the  second-law  analysis  in  a  rectangular 
duct  for  mixed  convection  regime.  The  simulation  of  mixed  con¬ 
vection  in  a  square  cavity  with  protruding  body  having  different 
aspect  ratios  was  performed  by  Yilbas  et  al.  [45].  Abu-Hijleh  and 
Heilen  [46]  studied  the  entropy  generation  due  to  laminar  mixed 
convection  from  an  isothermal  rotating  cylinder  for  three  cylinder 
radii  and  covered  wide  ranges  of  Reynolds  number  and  buoyancy 
parameter.  They  observed  that  entropy  generation  increased  as 
the  Reynolds  number  and  buoyancy  parameter  increased  and  it 
decreased  as  the  cylinder  radius  increased.  For  the  same  com¬ 
bination  of  Reynolds  number  and  buoyancy  parameter,  entropy 
generation  was  mainly  due  to  thermal  effects  at  small  cylinder  radii 
and  due  to  viscous  effects  at  large  cylinder  radii.  Abbassi  et  al.  [47] 
analyzed  the  entropy  generation  in  Poiseuille-Benard  channel  flow. 
He  found  that  the  maximum  entropy  generation  is  localized  at  areas 
where  heat  exchanged  between  the  walls  and  the  flow  is  maximum. 
There  is  no  significant  entropy  generation  is  seen  in  the  main  flow. 
In  the  similar  manner,  Nourullahi  et  al.  [48]  interested  by  the  issue 
of  entropy  generation  and  Nusselt  number  in  Poiseuille-Benard 
channel  flow  are  analyzed  by  solving  numerically  Navier-Stokes 
and  energy  equations  with  the  use  of  the  classic  Boussinesq  in 
compressible  approximation.  Their  result  showed  that  the  Nus¬ 
selt  number  changes  very  slightly  and  it  is  almost  constant  for 
low  values  of  inclination  angle.  It  decreases  for  higher  values  of 
inclination  angle.  The  entropy  generation  due  to  heat  transfer  is 
localized  at  areas  where  heat  exchanged  between  the  walls  and 
the  flow  has  a  maximum  value,  while  the  entropy  generation  due 
to  fluid  friction  is  maximum  at  areas  where  the  velocity  gradients 
are  maximum  such  as  vortex  centers.  In  other  words,  the  inclina¬ 
tion  angle  is  a  good  parameter  to  control  entropy  generation  in  the 
presence  of  buoyancy  forces.  Mahmud  and  Fraser  [49]  focused  on 
magnetohydrodynamic  mixed  convection  through  a  vertical  chan¬ 
nel  packed  with  fluid  saturated  porous  substances  by  applying  first 
and  second  laws  of  thermodynamics  to  analyze  the  problem.  They 
solved  governing  equations  analytically  under  reseasonable  sim¬ 
plifications  for  different  parameters  such  as  Hartmann  number, 
Plank  number,  Richardson  number  and  group  parameter  (Br//7). 
They  found  that  the  entropy  generation  number  is  characterized 
by  a  concave  shaped  profile  and  is  symmetric  about  the  chan¬ 
nel  centerline  for  a  symmetrical  temperature  boundary  condition. 
Radiation  and  mixed  convection  parameters  have  a  more  domi¬ 
nating  influence  on  entropy  generation  rate  than  porous-magnet 
and  group  parameters.  Expressions  for  the  idle  points  of  entropy 
generation  are  derived  and  their  location(s)  are  determined.  Irre¬ 
versibility  due  to  the  heat  transfer  and/or  fluid  friction  becomes 
zero  at  these  idle  points.  Shohel  and  Roydon  [50]  did  a  similar  work 
on  mixed  convection-radiation  interaction  in  a  vertical  channel  by 


calculating  entropy  generation.  Another  application  has  been  per¬ 
formed  by  Cheng  et  al.  [51]  on  entropy  production  in  mixed 
convective  flow  in  a  vertical  channel  with  transverse  fin  arrays. 

6.  Entropy  generation  of  natural  convection  in  nanofluid 
filled  enclosures 

Addition  of  nanoparticles  (A1203,  Cu,  Ag,  Au,  CuO,  Ti02,  etc.)  of 
any  material  into  base  fluid  such  as  water,  glycerol,  etc.  is  called 
as  nanofluid.  In  this  way,  heat  transfer  of  fluid  is  enhanced  due  to 
higher  thermal  conductivity  values  of  nanoparticles.  In  the  same 
time,  viscosity  of  fluid  is  changing  and  it  affects  the  heat  trans¬ 
fer  enhancement.  Studies  on  entropy  generation  of  heat  transfer 
in  nanofluid  filled  energy  systems  are  extremely  limited.  In  this 
context,  Shahi  et  al.  [65]  performed  a  numerical  study  on  entropy 
generation  due  to  natural  convection  of  a  nanofluid  that  consists 
of  water  and  Cu  in  a  cavity  with  a  protruded  heat  source.  They 
used  finite  volume  method  and  solved  the  equations  for  different 
Rayleigh  numbers,  solid  concentration  and  heat  source  location. 
Their  results  showed  that  the  maximum  value  of  Nusselt  number 
and  minimum  entropy  generation  are  obtained  when  heat  source 
mountains  in  the  bottom  horizontal  wall.  Singh  et  al.  [66]  present  a 
theoretical  investigation  of  the  entropy  generation  analysis  due  to 
flow  and  heat  transfer  in  nanofluids.  They  considered  the  most  com¬ 
mon  alumina-water  nanofluids  as  the  model  fluid.  They  observed 
that  at  lower  tube  diameter,  flow  friction  irreversibility  is  more  sig¬ 
nificant  and  at  higher  tube  diameter  thermal  irreversibility  is  more. 
For  both  laminar  and  turbulent  flow,  there  is  an  optimum  diameter 
at  which  the  entropy  generation  rate  is  the  minimum  for  a  given 
nanofluid. 

7.  Conclusions 

The  present  review  is  a  comprehensive  outlook  on  the  research 
progress  made  on  entropy  generation  in  natural  convection  and 
mixed  convection  in  ducts  or  channels.  The  important  result  can 
be  drawn  from  the  reviewed  literature  is  that  the  main  reasons  of 
entropy  generation  are  heat  transfer  and  fluid  friction  in  a  thermal 
system  with  buoyancy  force.  From  the  observed  results  it  is  clearly 
seen,  that  analysis  of  entropy  generation  is  crucial  for  all  thermal 
systems  to  reduce  energy  and  control  of  heat  transfer  and  fluid 
flow.  The  most  important  parameter  is  Rayleigh  number  on  entropy 
generation  in  natural  convection.  The  entropy  generation  increases 
with  increasing  of  Rayleigh  number  due  to  increasing  of  heat  trans¬ 
fer  and  fluid  friction.  Entropy  analysis  indicates  a  local  address 
that  where  energy  losses  high  in  the  physical  model  or  systems. 
Darcy  number  is  also  another  effective  parameter  for  entropy  gen¬ 
eration  in  porous  media  filled  systems.  The  entropy  also  increases 
with  increasing  of  Darcy  number.  In  the  case  of  mixed  convection, 
the  entropy  generation  decreases  with  increasing  of  Richardson 
number  due  to  decreasing  of  heat  transfer.  Hartmann  number 
is  another  effective  parameter  in  magnetohydrodynamics  studies 
with  electrically  conductive  fluids.  In  this  case,  increasing  of  Hart¬ 
mann  number  decreases  the  flow  velocity.  However,  the  entropy  is 
reduced.  Application  direction  of  magnetic  force  is  also  important 
to  minimize  the  entropy  generation.  Entropy  generation  decreases 
with  decreasing  of  irreversibility  ratio.  Inclination  angle  makes 
important  effect  in  natural  convection  studies  due  to  heat  transport 
way  changes  by  gravity.  In  these  kinds  of  studies,  entropy  genera¬ 
tion  is  also  affected.  The  geometrical  shape  and  aspect  ratio  of  the 
cavity  change  the  entropy  generation  even  at  constant  Rayleigh 
number.  And  type  of  the  fluid  with  Prandtl  number  places  impor¬ 
tant  role.  Recently,  the  second  law  thermodynamics  was  applied 
to  industrial  problems  to  minimize  the  entropy  generation  [60,61  ]. 
There  are  very  few  correlations  between  effective  parameters  and 
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entropy  generation.  Researcher  can  focus  on  that  subject  for  future 
studies.  The  trend  in  the  literature  goes  to  extension  of  natural  and 
mixed  convection  studies  by  adding  entropy  generation. 
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